LBNL-49070 



Statistical Errors in the Measurement of Particle Thresholds 

R. N. Cahn|] 

Lawrence Berkeley National Laboratory 
1 Cyclotron Rd., Berkeley, CA 94720 
(Dated: February 7, 2008) 

Simple rules of thumb are derived for the precision with which s-wave and p-wave thresholds can 
be determined by a series of equally spaced cross section measurements near threshold. Backgrounds 
and beam spread are ignored. 
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An important question for a linear collider is the precision with which particle thresholds can be measured. 
T-H , The thresholds in question may have either s-wave or p-wave behavior, that is, the cross section may rise as 
■ or p^. Measurements of the production cross section near the threshold determine particle masses quite 

' well, but the question is how well. The result must depend on the number of events registered, TV, and the 
interval over which the data are taken. The beamstrahlung will degrade the nearly monochromatic beam, but 
the effective beam energy spectrum still has a delta function at the nominal energy. The 50% or so of the beam 
that has has its energy reduced by 5% or so simply becomes a source of a diffuse background. Because the 
highest energy component of the beam occurs in the delta function, the sharp threshold in the s-wave survives. 
On For our purposes we ignore the beamstrahlung altogether. This means that our estimates are too optimistic. 
We also ignore the effect of the possible widths of the produced particles. Detailed threshold measurements 
might even measure these widths |2]. 
^ ' Here we suppose the energy is set at evenly spaced discrete values and a fixed amount of luminosity is 
I delivered. Let us suppose that at the ith energy value the expected number of events given the cross section 
— . with some trial values for the parameters and the delivered luminosity is n^^^. If the measured number of events 
' is 71°*"*, then probability of this set of measurements with these parameters is 
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o: ^-ne-<^^^ (1) 
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^ ' Let us suppose that we accumulate an integrated luminosity L at each of TV -f 1 energies, Eq -t- {i/N)AE, i = 

Qni 0, 1, ...TV. We suppose that the threshold occurs between Eq and Ei. We can take the log-likelihood function 

^ • £ = InP -I- constant to be a function of the number of observed events, n°^^, and the number of expected 
events, n^^^, at each energy: 



>< : C = f2 ^ Hnr/nf') - nf" + nf^ (2) 

S ; 

where we have added the n"'"' term so that C vanishes when n"'"' = n^^^. The values of n^^^ depend on the 
parameters describing the cross section. 

Now the expected uncertainty in a parameter is given by 
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evaluated at the true value of a. 

Here we shall be interested solely in the threshold for the process. The cross section can be written 

a = aa{{E~2m)lEthrY) (4) 

where 2m = Ethr is the trial value of the threshold. However, in order not to mix the unknown value of Ethr 
with the parameter m, we replace Ethr in the denominator by E'oi where Eq is the lowest energy setting, the 
one just below threshold. Thus Eq differs from Ethr typically by less than a percent. 
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We write 2m = Eq + \AE where m is the true mass and where < A < (1/A^). By hypothesis, the first 
non-zero contribution is for i = I. 
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where e = AE/Eq. 

In the p-wave case 1^ — 3/2 and we have 
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where we noted that doe^/^ = a{AE) is the cross section at the upper end of the AE interval and where 
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See Fig.0. Thus we find 
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Martyn and Blair ^ give an example of e^e^ -> [ibPr, with AE = 10 GeV, C7{AE) = 10 fb, TV = 10, L = 10 
fb~^. With these values we find 



a™ = 0.0745 GeY/ ^fp{\N,N) 



For A = 0. this gives 0.084 GeV, while Martyn and Blair give 0.09 GeV. 
In the s-wave case 
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We define 
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where the zeta function is defined by = S^i n so that /s(0, oo) = l See Fig. 0. We then find 

^™ =^^'/'^'/'C(3/2)/.(AiV,iV) 
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FIG. 1: The functions fp{z,N) and fs{z,N) for two values of A'^ 
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0.030 
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0.029 
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1.59 
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TABLE I: Comparison of Eq.(|l3|) and toy Monte Carlo study with A'^ = 10 data points for an s-wave threshold near 255 
GeV, with ao = 800 fb and 10 fb~^ at each point, corresponding to the study of Martyn and Blair for ej^e'^ XiXi ■ 
Martyn and Blair give an uncertainty of 0.04 GeV for a threshold of 255.0 GeV. As the threshold is placed closer and 
closer to the first data point, 256 GeV, the resolution steadily improves. The conservative estimate is obtained by taking 
A7V = 0. 
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where we noted that croe^^'^ = a-{AE) is the cross section at the upper end of the AE interval. 
For the example given by Martyn and Blair 0, i = 10 fb"\ a{AE) = 160 fb, = 10, we find 

a = 0.0276[/,(AA^, N)]-'^/^ (14) 



A comparison of a toy Monte Carlo study and Eq.(|l3|) is shown in the Table. 

The functions fs and fp reflect the possible fortuitousness of our scan. If we get a cross section measurement 
just above threshold, we get a much smaller error. The conservative errors are obtained by setting A = in 
these relations. We see that since the total luminosity used in the measurement is NL, the p-wave measurement 
doesn't depend directly on the choice of N, while the s-wave measurement benefits weakly from increasing N. 

With the approximations 
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fp{0,N) = 1-^ (15) 

we have the conservative estimates 
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